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It is argued that the observed vacuum energy density and the small values of the neutrino masses 
could be due to anthropic selection effects. Until now, these two quantities have been treated 
separately from each other and, in particular, anthropic predictions for the vacuum energy were 
made under the assumption of zero neutrino masses. Here we consider two cases. In the first, we 
calculate predictions for the vacuum energy for a fixed (generally non-zero) value of the neutrino 
mass. In the second we allow both quantities to vary from one part of the universe to another. 
We find that the anthropic predictions for the vacuum energy density are in a better agreement 
with observations when one allows for non-zero neutrino masses. We also find that the individual 
distributions for the vacuum energy and the neutrino masses are reasonably robust and do not 
change drastically when one adds the other variable. 



I. INTRODUCTION 

The smallness of the vacuum energy density pa (the 
cosmological constant) is one of the most perplexing mys- 
teries that we now face in theoretical physics. The val- 
ues of pa indicated by particle physics models are far too 
large, at least by 60 orders of magnitude. The only plau- 
sible explanation that has so far been suggested is that 
PA is a stochastic variable, taking different values in dif- 
ferent parts of the Universe, and that the smallness of the 
observed value is due to anthropic selection effects. The 
proposed selection mechanism is very simple* [1—8]. The 
growth of density fluctuations leading to galaxy forma- 
tion effectively stops when the vacuum energy dominates 
the Universe. In regions where p\ is greater, it will dom- 
inate earlier, and thus there will be fewer galaxies (and 
therefore fewer observers). Suppose p\ is the value for 
which the vacuum energy dominates at about the epoch 
of galaxy formation, za ~ Zg- Then, values pa 3> /5a 
will be rarely observed, because the density of galaxies 
in the corresponding regions is very low. On the other 
hand, values p\ -C pa are also rather unlikely, simply 
because this range of values is very small. (This argu- 
ment assumes that the prior distribution for pa is not 
sharply peaked at pa = 0; see below.) Thus, a typical 
observer should expect to measure pa ~ Pa- the vacuum 
energy should dominate roughly at the epoch of galaxy 
formation. 

This anthropic prediction was subsequently confirmed 
by supernovae and CMB observations. The observed 
value [12,13] for the vacuum energy density, in units of 



the critical density, fl A bs as 0.7, is in a qualitative agree- 
ment with the anthropic probability distribution, which 
is peaked at [2,3,6,7] J7a ~ 0.9. The most recent analysis 
[8] finds that anthropic predictions are in agreement with 
observations at a 2a level. 

Here we show that the agreement can be further im- 
proved by allowing for a nonzero neutrino mass. The 
value pa at the peak of the probability distribution for 
Pa depends on the magnitude of density fluctuations on 
the galactic scale, <jq. The calculation in [8] was per- 
formed assuming massless neutrinos and a scale- invariant 
spectrum of fluctuations, normalized using the large-scale 
CMB measurements. If uq is decreased, galaxies form 
later, and pa becomes smaller, thus improving the agree- 
ment with p°a s ■ This is precisely the effect of massive 
neutrinos, since they slow down the growth of density 
fluctuations. If the sum of the three neutrino masses is 
m v ^1 — 2 eV, this suppression is sufficiently large to 
bring the anthropic prediction for pa into agreement with 
observations at the la level. 

It is very intriguing that the value m v ~ 1 eV for 
the sum of the neutrino masses was itself recently pre- 
dicted from anthropic considerations. It was suggested 
in [14] that the neutrino masses might also be stochas- 
tic variables, changing from one part of the Universe to 
another, and that their small values might be due to an- 
thropic selection. In regions with a larger value of m„, 
the suppression of density fluctuations is stronger. As a 
result, there are fewer galaxies in such regions, and the 
corresponding values of m„ are less likely to be observed. 
The probability distribution for m u calculated in [14] is 
peaked at m„ as 2 eV, with a la range between 0.5 eV 
and 3 eV. This corresponds to a mild but non-negligible 



*For earlier suggestions that the cosmological constant may 
be anthropically selected, see [9-11]. 
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suppression of galaxy formation. ^ Indeed, recent analy- 
sis of X-ray cluster data combined with CMB and 2dF 
has suggested that m v = 0.64^g' 2 g eV [16], although this 
remains controversial given cluster-related systematic un- 
certainties [31]. In this paper, we shall also consider mod- 
els where both m u and p\ are assumed to be stochastic 
variables. 

The paper is organized as follows. In Section II we 
describe our general approach and justify our choice of 
the theoretical prior distributions for m v and p\. In Sec- 
tion III we calculate the probability distribution for pa 
at nonzero values of m„, as well as the joint probabil- 
ity distribution for m v and p\ in models where both are 
variable. We conclude with a discussion of our results 
in Section IV. Some technical details are given in the 
Appendix. 



II. THE METHOD AND THE PRIORS 

The method used in this paper follows closely that in 
[7,14]. Consider a model in which m v and pA are allowed 
to vary from one part of the Universe to another. We 
define the probability distribution V '{m v , p ;C)dm v dp a as 
being proportional to the number of observers in the Uni- 
verse who will measure m„ in the interval dm u and pa in 
the interval dp a- This distribution can be represented as 
[2,7,14] 

OC T^prior{jFiv •> Pa) 

x f dM n G (m v ,pK,M) N obs {m v ,pK,M). (1) 

Here j ^PpT', 

which is proportional to the comoving volume of those 
parts of the Universe where m„ and p\ take values in 
intervals dm v and dp\, n G (m Vl p\, M)dM is the comov- 
ing number density of galaxies of mass in the interval 
dM that will ever form in regions with given values of 
m„ and p\, and N b s (m u , p\, M) is the average number 
of observers per galaxy. The distribution (1) gives the 
probability that a randomly selected observer is located 
in a region where the sum of the three neutrino masses 



\m v , pAjdm^dpA is the prior distribution, 



tit is conceivable that the suppression of galactic-scale den- 
sity fluctuations relative to the largest observable scales is not 
due to massive neutrinos, but reflects some feature in the pri- 
mordial fluctuation spectrum. Such suppressed small-scale 
fluctuations have been suggested by certain galaxy cluster 
observations (assuming that the matter density parameter is 
!)jm ~ 0.3, see [15]), and by the WMAP team [13] who pointed 
out that an initial spectrum with a decreasing running spec- 
tral index is somewhat preferred by the data. Whatever the 
origin of such a suppression, it improves the agreement of the 
anthropic predictions for pA with observations. 



is in the interval dm v and the vacuum energy density is 
in the interval dp a- 

We shall assume, as it was done in [6,7,14], that the 
integral in (1) is dominated by large galaxies like the 
Milky Way, with mass M > M G ~ 10 12 M Q . We shall 
also assume, as a rough approximation, that for galaxies 
in this mass range N b a does not depend significantly on 
to„ or pl and is simply proportional to the number of 
stars, or to the mass of the galaxy, 



N obs (m v ,p A ,M) cx M . 



(2) 



(We shall comment on the validity of this assumption in 
Section IV.) Then it follows from Eq. (1) that 

V(m v ,p A ) oc V P rior(m v ,pA)F(M > M G ,m v ,p A ), (3) 

where F(M > M G , m„, pa) is the fraction of matter that 
clusters into objects of mass larger than M G in regions 
with these values of m„ and pa- 

The fraction of collapsed matter F(M > M G ,m u , pa) 
can be approximated using the standard Press-Schechtcr 
formalism [17]. We assume a Gaussian density fluctua- 
tion field (5(x, t, m u , pa) with a variance <r(t, m l/ ,pA) on 
the galactic scale (M G = 1O 12 M ), 



P{8,t,m l/ , pa) oc exp 



2c 2 (t,m I/ ,pA) 



(4) 



At early times, when neutrinos are still relativistic, the 
variance a is assumed to be independent of m v and pA- 
An overdense region collapses when the linearized den- 
sity contrast 5 exceeds a critical value <5 C , which in the 
spherical collapse model with massless neutrinos is esti- 
mated to be [1] 5 C ss 1.63. We shall assume that the same 
value applies when neutrinos have a small mass. 

In the Press-Schechter approximation, the asymptotic 
mass fraction in galaxies of mass M > M G is 



F(M > M G , m v , pa) oc crfc 



X 



V2 



where 



X = S c /<Joo(m u , pa) 
and ^(m,,, p A ) = a(t — > g»,to„,pa)- 



(5) 



(6) 



A. Prior distributions 

The calculation of the prior distribution P pr ior(m u , pl) 
requires a particle physics model which allows the neu- 
trino masses and the dark energy density to vary and 
a cosmological "multiverse" model that would generate 
an ensemble of sub-universes with different values of ra v 
and pa- Let us first focus on the neutrino masses. (The 
discussion in this subsection follows [18,14].) 
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Dirac-type neutrino masses can be generated if the 
Standard Model neutrinos v a mix through the Higgs dou- 
blet VEV $ to some gauge-singlet fermions u%, 

g afl *P a v§. (7) 

The couplings g a p will generally be variable in string the- 
ory inspired models involving antisymmetric form fields 
F a interacting with branes. (Here, the index a labels dif- 
ferent form fields.) F a changes its value by AF a = q a 
across a brane, where q a is the brane charge. In the low- 
energy effective theory, the Yukawa couplings g a p become 
functions of the form fields, 

g a p-2^9 a f3 th — + ••• ( 8 ) 

lVlp 

Here, the summation is over all form fields, the coeffi- 
cients <7®£ are assumed to be numbers ~ 1 , and M p is the 
effective cutoff scale, which we assume to be the Planck 
mass. We have assumed that g a p vanish at some point 

F a = Fa°^ in the f-space; then Eq. (8) is an expansion 
about that point. 

In such models, closed brane bubbles nucleate and ex- 
pand during inflation [19], creating exponentially large 
regions with different values of the neutrino masses. 
When F a changes in increments of q a , m v changes in in- 
crements of Am„ ~ Qq a /Mp. To be able to account for 
neutrino masses < 1 eV, we have to require that Ato„ < 1 
eV, that is, 

q a < lO~ n Ml (9) 

for at least some of the brane charges. Such small val- 
ues of the charges can be achieved using the mechanisms 
discussed in [20-22]. 

The natural range of variation of F a in Eq. (8) is the 
Planck scale, and the corresponding range of the neutrino 
masses is < mi'' < $. (Here, the index i labels the 
three neutrino mass matrix eigenvalues.) Only a small 
fraction of this range corresponds to values of anthropic 
interest, m v < 10 eV. In this narrow anthropic range, 
we expect that the probability distribution for F a after 
inflation is nearly flat [23], 

dVprior oc dF 1 dF 2 ..., (10) 

and that the functions g a p(F a ) are well approximated by 
linear functions (8). If all three neutrino masses vary 
independently this implies that 

dVpriorOtdm^dml^dmP. (11) 

The probability for the combined mass m v = mty to 
be between m v and m v + dm v is then proportional to the 
volume of the triangular slab of thickness <~ dm v in the 
3-dimensional mass space, 

dVp r i or oc m 2 ,dm l ,. (12) 



Alternatively, the neutrino masses can be related to one 
another, for example, by a spontaneously broken family 
symmetry. If all three masses are proportional to a single 
variable mass parameter, then we expect 

dVprior oc dm v . (13) 

To reiterate, the prior distributions (12), (13) arc 
generic if the following two conditions are satisfied: (i) 
the full range of variation of m v is much greater than the 
anthropic range of interest and (ii) the point m v = is 
not a singular point of the distribution. The values of F a 
which make the Yukawa couplings (8) vanish are not in 
any way special from the point of view of the distribu- 
tion for F a , so there is no reason to expect a singularity 
at that point, and thus the condition (ii) is likely to be 
satisfied. 

As was shown in [14], the anthropic prediction for the 
neutrino masses with the prior given by Eq. (12) is in a 
strong disagreement with current observational bounds. 
Allowing for a varying pa is unlikely to change this result. 

There is also a possibility that the right-handed neu- 
trinos i/P in (7) have a large Majorana mass Mr » In 
this case, small neutrino masses can be generated through 
the see-saw mechanism, 

m v ~ g 2 <S> 2 /M R . (14) 

If Mr is variable, say, within a range Mr < M p , then its 
most probable values are likely to be ~ M p , and the prior 
distribution will favor m v ~ 10" 6 eV. Such extremely 
small values of m v are in conflict with the neutrino oscil- 
lation data. In this paper we shall assume Dirac masses 
with the prior (13). 

It should be noted that the Higgs potential and the 
Higgs expectation value $ in (7) will generally be func- 
tions of F a . Moreover, each field F a contributes a term 
F 2 /2 to the vacuum energy density pa, and regions with 
different values of F a will generally have different values 
of pa- However, in the presence of several form fields 
with sufficiently small charges, variations of all these pa- 
rameters are not necessarily correlated, and here we shall 
assume that there is enough form fields to allow indepen- 
dent variation of the relevant parameters. We can then 
consider sub-ensembles of regions where some of the pa- 
rameters are variable, while the other are fixed. 

In this paper we will be concerned with two models: (i) 
variable pa at a fixed value of m v and (ii) variable m v and 
Pa, all other parameters being fixed. A line of reasoning 
similar to that above suggests that, in the anthropically 
interesting range, the prior distribution for pa should also 
be flat [4,5] 

dV prior oc dp a- (15) 

Assuming independent variation of m„ and pa, we then 
have 

Vp r i or {m v ,pA) = const. (16) 



3 



We finally comment on the recent work [24-26] sug- 
gesting that string theory involves a large number of form 
fields F a (~ 10 3 ) and therefore admits an incredibly large 
number of vacua (~ 10 1000 ) characterized by different 
particle physics parameters.-'' The spectrum of p\ and 
m v could then be very dense even if the brane charges 
are not small, q a <~ M^. However, in this "discretuum" 
of vacua, nearby values of m v or pa correspond to very 
different values of the form fields, and we can no longer 
argue that the prior distribution should be flat in the an- 
thropic range. Calculation of the prior in such a discre- 
tuum remains an important problem for future research. 
In this paper we shall use the distribution (16). 



B. The full distribution 

The full probability distribution can now be obtained 
from Eqs. (3), (5), and (16), 



where 



V(m v ,pA) oc erfc 
where X is given by Eq. (6), 

X = Sc/doo- 



X 



V2. 



(17) 



(18) 



Depending on the model, the asymptotic density contrast 
(Too is a function of m v and pa or of ra v only. 



III. THE PROBABILITY DISTRIBUTION 

Let us consider models where both p\ and m v are al- 
lowed to vary from one part of the Universe to another. 
We shall calculate the probability distribution expressed 
in terms of p\. We will also state some of our results in 
terms of 



n A = {8irG/3)p A h- 2 = uj A h- 



(19) 



To describe the growth of density fluctuations, it will 
be convenient to use the ratio 



PA 



PAl(t) 



(1 + *)- 



(20) 



as our time variable. The variance c7 co (m I/ ,pA) in 
Eq. (18) is proportional to the total linear growth fac- 
tor of the density perturbations, 



D°°(U,pa) = D(f v ,p A ,x^> oo) , 



(21) 



"'"Apart from the form fields, the string vacua are character- 
ized by a number of other parameters specifying the string 
compactifications. 



dm 



(22) 



is the neutrino fraction of the total matter density. 
Hence, the effect of neutrino masses on the fraction 
F(M > MG,m„ 7 pa) in (17) can be taken into account 
by writing 



^ooK.Pa) = CToo(0,/0a) 



£°°(0,pa) 



(23) 



where D°°(0,pa) is the total growth factor in regions 
with /„ = 0. 

The linear growth factor as a function of redshift can 
be found numerically, using standard codes such as CMB- 
FAST [27], or analytically, using, e. g. fitting formulae 
from [28] . We have tested that the two methods give re- 
sults that agree within a margin of error of a few percent. 
Technical details related to their evaluation can be found 
in the Appendix. Throughout this paper we choose to 
work with the semi-analytical fitting formulae given in 
Sec. A 2. 

We need to relate the asymptotic density contrast cri? 
in Eq. (23) to the locally observable quantities. Current 
CMB measurements do not extend down to the galactic 
scale, and to sidestep issues related to light-to-mass bias, 
the galactic density fluctuations are usually inferred by 
extrapolating from larger scales (on which the effect of 
neutrino free streaming is unimportant) using a nearly 
scale-invariant fluctuation spectrum and assuming that 
the neutrino masses are negligible. The growth factor 
evaluated in this way overlooks the non-zero massive neu- 
trino fraction /* in our part of the Universe, which is cur- 
rently unknown but potentially significant. This growth 
factor, D(f*,p A ,x), is evaluated under the assumption 
that neutrinos are massless, but with f2* added to £l c dm- 
Here and below, hats indicate the quantities evaluated for 
our local region assuming massless neutrinos, and super- 
script "*" indicates the actual local values of these quan- 
tities. For example, x* w 2.3 is the present local value of 
x and a is the asymptotic density contrast in our region 
calculated assuming massless neutrinos. Hence, we can 
write 



Ooo (0) = &c 



£°°(0,Pa) 



(24) 



where D°°(K,p\) = D(f:,p* A ,x^ oo). 

To evaluate (Too, we again note that in the absence 
of massive neutrinos, the growing mode of density fluc- 
tuations is proportional to D(0, p A ,x). Therefore, the 
asymptotic density contrast aoc is related to the present 
value a(M) via 



(Toe = <t(M) 



£°°(/;,Pa) 
D(f*,p* A ,x*) 



(25) 
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Finally, using Eqs. (23), (24) and (25), we can rewrite 
Eq. (18) as 



X 



5 C D(f;,p* A ,x*) 
a(M) £><"(/„, p A ) 



(26) 



Here, <r(M) is the density contrast on the galactic scale 
M = 10 12 M© inferred from the large-scale CMB data as- 
suming massless neutrinos. We emphasize that the actual 
density contrast should be smaller if indeed the suppres- 
sion of the galactic density fluctuations due to neutrinos 
is non-negligible.'*. To estimate a{M) we can use the lin- 
ear power spectrum inferred from measurements of CMB 
anisotropics. First we need to find the length scale Rq 
corresponding to the mass scale M = 10 12 M Q . It can be 
found from 



R(M) = 



( 3M 

V47T/90 



1/3 



0.951 Q]^- 3 h- 2/3 



M 



10 12 M m 



1/3 



Mpc , (27) 



where po » 1.88 • 10 26 il^jh 2 kg/m 3 is the mean cosmic 
density. For h — 0.72 and VLm = 0.27 this gives 



R G « 1.32 h _1 Mpc . 



(28) 



The corresponding linearized density contrast gr g found 
using WMAP's best fit power law model [13] is 

gr g w<r(10 12 Af Q ) w 2.41 ± 0.26 . (29) 

Substituting analytical expressions for D(f*, p* A ,x*) 
and D°°{f Ul pa), derived in the Appendix (Sec. A 2), into 
Eq. (26) and using the fact that x* » x eq gives 



X 



^ ( * \- 

1-/^ 4P/3 

i -ft 



1/3 



1-p 



p/3 



[G(oo)]p 



where 



V25 - 24/, - 1 



(30) 



(31) 



and G(x) is the growth factor in a flat universe filled with 
pressureless matter and vacuum energy, given by [29,6]: 



G(x) 



5 


"1 +x~ 


1/2 r dw 


6 


X 


J W 1 /6(l + M; )3/2 



(32) 



^When better data for the actual linearized galactic-scale 
density contrast is available, it can also be used for the evalu- 
ation of F(M > Mo,rn v , pa)- The only modifications needed 
in Eq. (26) are that a(M) has to be replaced by a*(M) and 
the factor D{f* , p\,x*) should be replaced by D(f*,p A ,x*). 



with the asymptotic value G(oo) w 1.44. 

The probability distribution V(m u ,p\) is obtained 
from Eqs. (17) and (26) and depends on parameters p* , 
(7fl G and /*. The last parameter, the value of /„ in our 
part of the Universe, is one of the quantities we are trying 
to predict. On the other hand, we found that the distri- 
bution is not very sensitive to the assumed value m*: 
varying m* between 0.1 and 3 eV changes the location of 
the peak of distribution functions by no more than 10% 
percent. This uncertainty can be accounted for by multi- 
plying P(m„, pa) by an appropriately chosen prior proba- 
bility p(m*), reflecting the current state of our knowledge 
of the value of m*, and integrating over m*. Current 95% 
upper bound on the sum of neutrino masses, obtained 
using the WMAP and SDSS data [31], is m* < 1.7 eV. 
Hence, we choose p(m*) to be a Gaussian centered at 
to* = with a standard deviation of o m > = 0.85eU. A 
similar procedure is also needed to account for the un- 
certainty in our knowledge of p\ and ctr g . We compute 



fl* A h* 



using h* 



h±a h = 0.72 ± 0.08 measured 
by the HST project [30] and tt* A = tt A ±<rn A = 0.7±0.08, 
a conservative estimate from CMB and large-scale struc- 
ture observations [13,31] . We take p(w A ) to be a Gaus- 
sian centered at = ^A" , with a standard deviation 
a UJA = Cl\2hah + han A . We account for the uncertainty 
in <tr g taking p(aR G ) to be a Gaussian with the mean 
and the standard deviation given by Eq. (29). Hence, we 
have 



V{m v ,p A )^ J dm*p(m*) J dw* A p(uj* A ) 

x / da RG p(a Ra ) V(m u ,pA) . (33) 




0.01 0.1 



1 10 
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100 1000 
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FIG. 1. The logarithmic probability distributions pjyV(m„, pa) 
vs Pa/Pa f° r three fixed values of m v : 0, 1 and 2 eV, peaked from 
right to left. In each of the three cases, the lightly and densely 
shaded areas are the regions excluded at 68% and 95% level respec- 
tively. Also shown (the dark vertical strip), is the la uncertainty 
in p A derived from observations. 

In the following two subsections we consider two cases. 
In the first, we calculate the probability density for p\, 
while using several fixed values of m v . In the second, 
we let both p\ and m v vary and calculate their joint 
probability as well as the effective probability for each of 
the two parameters. 



A. P{pa) for fixed values of m v 

In this subsection we consider a few fixed values of m„ 
and evaluate corresponding probabilities for pa- 

The probability P(pa) per log p\ for m„ = 0, 1 and 
2 eV is shown in Fig. 1 along with the 1 and 2a bounds 
in the three cases. These bounds are found by cutting 
off equal measures of probability on both ends (16% and 
2.5% for 1 and 2a, respectively). Also shown in the figure 
is the observationally favored range of p\. As one can see 
from the figure, the observed value is within the 2a range 
of the anthropic prediction in all three cases and within 
the la bound in the m v = 2eV case. 



B. Joint probability P(m v ,pA) 

Let us now calculate the probability distribution 
V{m v ,p\) with both m v and pa variable. 




We are interested in finding the 68% and 95% probabil- 
ity contours in the m„-pA parameter plane. The simple 
procedure of cutting off equal measures of probability 
on both ends used in the previous subsection does not 
readily generalize to the case of more than one variable. 
Instead, we find contours of constant logarithmic proba- 
bility density m v p\P{m u , p\) that bound 68% and 95% 
of the volume under the surface defined by the probabil- 
ity distribution "P(m„, /ja). These contours are shown in 
Fig. 2. 



E 

— o.: 




FIG. 3. Effective logarithmic probability distribution for 
ra v and the boundaries of the 68% and 95% probability re- 
gions. 

The ensemble averaged values of the dark energy den- 
sity and the neutrino masses found using the joint prob- 
ability distribution V(m u , pa) are 

(fi A ) = 0.85 , {pa/ Pa) = 21 . M = l-5eV . (34) 



FIG. 2. The 68% and 95% probability contours in the m v -p^ 
parameter plane for flat priors. The shaded vertical strip shows the 
IcT uncertainty in the observed value of pjy. 
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FIG. 4. The effective probability from Eq. (36) plotted per 
log pa/ Pa- The lightly and densly shaded regions are excluded 
at 1 and 2a level, respectively. Expressed in terms of £2a, the 
la region is given by 0.695 < fi A < 0.988. 

Next we find the effective one dimensional probability 
distributions for m„ and p\. The probability distribution 
for m„ is found by integrating V(p\, m v ) over p\: 

V e s{m v )= I dptJ>{pK,m v ) . (35) 

In Fig. 3 we show the plot of V e s(m v ) per log m v , along 
with the boundaries of the 68% and 95% probability re- 
gions. 

In Fig. 4 we show 

V e s(pA) = ^ dm v V{pK,m v ) , (36) 

plotted per log p\, with the 68% and 95% probability 
regions. We find, as might be expected, that allowing for 
variable neutrino masses has a significant effect on the 
anthropic prediction for p\. In particular, it results in 
the currently observed value, f^A = 0.7 ±0.08, now being 
within the ler allowed region. 

We note that the integrated distributions (35), (36) 
should be interpreted with care. In Eq. (35), for exam- 
ple, the integration is performed over the full range of pA, 
which includes a large range of values which are already 
observationally excluded. Hence, this distribution can- 
not be used directly to make observational predictions 
for m v in our region. 

IV. DISCUSSION 

In conclusion, we have argued that both the small 
vacuum energy density and small values of the neutrino 



masses may be due to anthropic selection. We have con- 
sidered a model where m v and p\ vary independently, 
with flat priors for both variables. The predictions of 
this model are qualitatively similar to those derived by 
considering sub-ensembles with only one of these param- 
eters variable and the other fixed. But there is an in- 
teresting quantitative difference. The probability distri- 
bution for the joint ensemble is shifted towards smaller 
values of both m„ and p\. As a consequence, the agree- 
ment between the predicted and observed values of p\ is 
improved. 

In this work we have only considered the case of posi- 
tive p\. Negative values of p\ were discussed in [7] and in 
[32]. There, it was found that the probability of p\ < 
is comparable to that of p\ > 0. 

In our analysis we have assumed that the number of 
observers does not sensitively depend of m v or /?a and is 
simply proportional to the mass of the galaxy, Eq. (2). 
This appears to be a reasonable assumption for small 
values of m v and p\, such that they have little effect 
on structure formation. For larger values, the properties 
of the galaxies, including the number of habitable stars 
per unit mass, may be affected. For example, if p\ gets 
large, structure formation ends early and galaxies have a 
higher density of matter. This may increase the danger 
of nearby supernova explosions and the rate of near en- 
counters with stars, large molecular clouds, or dark mat- 
ter clumps. Gravitational perturbations of planetary sys- 
tems in such encounters could send a rain of comets from 
the Oort-type cloud towards the inner planets, causing 
mass extinctions. Further increase of density may lead to 
frequent disruptions of planetary orbits [33]. The effect 
of all this is to suppress the probability density at large 
values of p\ and thus to shift the distribution towards 
smaller values. This may further improve the agreement 
of the anthropic predictions with the data. With better 
understanding of galaxy and star formation, we may be 
able to obtain a quantitative measure of these effects in 
not so distant future. 
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APPENDIX A 

To compare the growth of linear density fluctuations 
in models with different values of m v and p^ 1 we choose 
the parameters so that all models have the same density 
contrast evolution while neutrinos arc still rclativistic and 
dark energy is negligible. This is achieved by fixing the 
values of w c = il c dmh 2 7 uib = &bh 2 and Tqmb today. 
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Then, for given values of /„ and pa = 3uja/(8ttG), we 
find 



UJ c +U) b 

and obtain fife, C( i m , J7a and f2„ using 
o-^io - ^£ o - 



(Al) 



fv{^cdm + ^b) 



(A2) 



We use the values a;;, = 0.024 and u c = 0.132, estimated 
from CMB and large-scale structure observations [13,31]. 



1. Evaluating growth functions using CMBFAST 

The matter growth functions can be evaluated nu- 
merically, using CMBFAST [27]. Given the values of 
cosmological parameters, including /„ and p\, CMB- 
FAST can compute Fourier components of the matter 
density contrast, 5{f Ul pA,k,x), where x is defined in 
Eq. (20). On scales much smaller than the neutrino 
free-streaming scale, k 3> kf a , neutrinos do not partic- 
ipate in gravitational clustering, so one expects ratios 
such as S(fl, p' A ,k,x')/5(fv, p\,k,x) to be independent 
of k. Galactic scales, corresponding to k g ~ IMpc -1 , 
satisfy this requirement for m v <C 30 eV, which includes 
the range of interest to us **. Therefore, since Eq. (26) 
contains ratios of growth factors, we can simply replace 
D(f*,p A *,x*) and D°°(f v ,p A ) by S(f*,p A *,k g ,x*) and 
$(fi>i PA,k g ,x — > oo). Namely, we write 



X = 



&(M) 5{j v , pa, k g ,x — > oo) 



(A3) 



where 5(f*,p\*, k gi x*) is the density contrast evaluated 
assuming neutrinos are massless, but with Q,* added to 
O c rf m . We have checked that, as expected, the results 
are not sensitive to a particular choice of the value of k, 
provided that k is sufficiently large. 



2. Evaluating growth functions analytically 

In addition to numerical methods, there are analyti- 
cal fitting formulae describing the evolution of density 
fluctuations in the presence of massive neutrinos. The 
growth of linear density fluctuations on scales k » kf s in 



*Tn addition to this requirement, in order for this ratio of 8s 
to be independent of k, the scale under consideration had to 
enter the horizon sufficiently long before z eq . Galactic scales, 
corresponding to k g ~ IMpc" 1 satisfy this constraint as well. 



a Universe dominated by nonrelativistic matter is given 
by [34] 



<t(z) cx (1 + z)-p, 



where 



p=\ (V25 - 24/„ - l) « 1 - jj/„, 



(A4) 



(A5) 



and the last step in (A5) assumes that f v <C 1. The 
growth of fluctuations effectively begins at the time of 
matter domination and terminates at vacuum domina- 
tion. 

The effect of massive neutrinos on the growth factor 
can be calculated using [28] 



D{f v , P A,x) = [D(f v ,p A ,x)]P. 



(A6) 



The functional dependence on /„ in Eq. (A6) is not ac- 
curate for very small values of m„, when neutrinos be- 
come nonrelativistic well after matter-radiation equality. 
But in this case the /^-dependence is very weak anyway. 
It has been verified in [28] that Eq.(A6) is accurate to 
within 3% in the whole relevant range of neutrino masses. 

In a flat universe filled with pressureless matter and 
vacuum energy the growth factor, G(x), is given by 
Eq. (32). For x <C 1, it reduces to the familiar growth 
function in a matter-dominated universe, G(x) w x 1 / 3 cx 
(l + z)-\ 

The effect of radiation can be included using the ex- 
act formula for a matter plus radiation universe (in the 
absence of dark energy) [35], 



G(x) = l + ±(x/x eq ) 1 ' 3 



(A7) 



The density of radiation is negligible when the dark en- 
ergy becomes important, and vice versa. Hence, we can 
write, to a good accuracy, 



D(.U,x) = l + -x- q ^(.U)G(x) , 



where 



x eq ^ {fv) — ( x eq) 



-1/3 



1 - ft 



4/3 



UIA 



1/3 



(A8) 



(A9) 



Note that the growth factor is normalized so that 
D(f v ,x — » 0) = 1. Equation (A9) follows from 

z" 1/3 = (1 + z eq )/(l + za) cx <Vfi- 1/3 

oc (1 - /,r 4/3 K + ucf'^l 1 ' 3 > (A10) 
and, since uib and uj c are fixed, we have 

Z-V3 K(1 _ /w) -4/3 w -l/3. (An) 

Eq. (A9) also depends on x* , currently estimated to be 
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-1/3 



2820 



(A12) 



and /*, the value of which is currently unknown. As 
discussed in Sec. Ill, the dependence of our results on /* 
is rather weak and we opt to marginalize over it. 

An approximation for G(x) accurate for all values of 
(x) was derived in [14]: 



G{x) 



,1/3 



1 



G(oof 



-\ -l/3c 



(A13) 



where a = 0.795. At large x, the growth of density fluc- 
tuations is stalled, and G(x) approaches the asymptotic 
value G(oo) w 1.44. 
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